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The asymmetric resonant exchange qubit under the influence of electrical noise 
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We investigate the influence of electrical charge noise on a resonant exchange (RX) qubit in a 
triple quantum dot. This RX qubit is a variation of the exchange-only spin qubit which responds 
to a narrow-band resonant frequency. Our noise model includes uncorrelated charge noise in each 
quantum dot giving rise to two independent (noisy) bias parameters £ and A. We calculate the 
energy splitting of the two qubit states as a function of these two bias detuning parameters to find 
“sweet spots”, where the qubit is least susceptible to noise. Our investigation shows that such 
sweet spots exist within the low bias regime, in which the bias detuning parameters have the same 
magnitude as the hopping parameters. The location of the sweet spots in the (e. A) plane depends 
on the hopping strength and asymmetry between the quantum dots. In the regime of weak charge 
noise, we identify a new favorable operating regime for the RX qubit based on these sweet spots. 


I. INTRODUCTION 

Universal cmantum computing with electron spins in 
quantum dotflhas been investigated within a wide range 
of possible implementations in solid state systems. Many 
implementations consider quantum dots in semiconduc¬ 
tors, mostly GaA^^ or silicoiP. One common feature 
of these implementations is their dependence on the con¬ 
trol of electric and magnetic fields on the nanometer scale 
which is needed for universal spin control^. This depen¬ 
dence couples the qubits to electric or magnetic nois^ 
introduced by the electric field of the gate voltages, the 
externally applied magnetic field, and the fluctuating ef¬ 
fective magnetic field produced by the nuclear spins. Al¬ 
though spin coherence times can be fairly long, the sus¬ 
ceptibility to electromagnetic noise implies a limitation 
of the coherence tim^. Therefore, we are interested in 
an implementation of a qubit which is protected against 
noise, in addition to faster qubit control techniques to 
enable as many coherent operations as possible within 
the spin lifetime. 

Qubits realized with the spin 1/2 of a single electron 
can be controlled using pulsed magnetic or (via spin-orbit 
coupling) electric oscillatory fields, and the exchange in¬ 
teraction can be used to couple such qubits to perform 
two-qubit gate^Il^. Encoded qubits can be operated with 
a reduced amount of magnetic or spin-orbit control. The 
encoding into the spin singlet and one of the spin triplet 
states of two electrons each localized on one of two nearby 
quantum dots allows for single qubit rotations generated 
by the exchange interaction and a static magnetic field 
gradienfPHni. Readout and spin preparation is imple¬ 
mented in singlet-triplet qubits by lowering the potential 
of one dot and transforming the spin information effec¬ 
tively into a charge signal by means of the Pauli princi¬ 
ple. Tw o-qubit o perations are possible with electrostatic 
or exchanges! between quantum dots be¬ 
longing to different qubits. 

Ultimately, a three electron-spin encoding allows for 
full qubit control without any magnetic fields and with¬ 
out relying on the spin-orbit coupling, but with the elec¬ 
trically controllable exchange interaction onl}^^^^^. The 


exchange-only scheme requires the exchange coupling be¬ 
tween pairs of spins to be switched on only for a short 
period of time. The fact that the exchange coupling 
can be switched off whenever the qubit is idle allows 
for an advantageos isolation of the spin qubit from the 
surrounding charge noise. The exchange-only qubit can 
be supplemented with additional control using an oscilla¬ 
tory (typically radio-frequency) electric fie ld wh en the ex¬ 
change interaction is constantly turned on.However, 
this enhanced control comes with additional decoherence 
channels due to the effects of charge noise, because the 
spin singlet and triplet states have slightly different or¬ 
bital wavefunctions when the exchange coupling is turned 
on. A substantial amount of experimental research has 
been done on the exchange-only spin qubit implemen¬ 
tation since its discovery, e.g. coherent control of the 
qubilPIl and readout and qubit preparatiorP^. For read¬ 
out and spin preparation, techniques used previous ly for 
singlet-triplet double quantum dots can be adaptecP^^. 
The decoherence produced by the hyperfine interaction 
with a nuclear bath and by electron-phonon interactions 
turn out to be of similar magnitude for the subspace 
and subsystem qubit^^. Here, we consider the influ¬ 
ence of charge noise, e.g. from the gate electrodeJ^ or 
from impurities in the materiaP^. It has been predicted 
that charg e noise can indeed be the dominant source of 
noise 

A strength of the exchange qubit with the exchange 
interactions switched on permanently is the suppression 
of low-frequency noise, giving rise to a regime in which 
the system responds to a resonant, narrow frequency 
bancpSHHMl^ This so called resonant exchange (RX) 
qubit allows for single-qubit control performed by ra¬ 
dio frequency signals in the resonant freque ncy in stead 
of pulse sequences of the exchange interact iorP^^^. Due 
to this narrow-band response, one can expect a natu¬ 
ral protection of the qubit against low frequency elec¬ 
tric charge noise. This expectation was confirmed in the 
case of one-parameter electric charge noise in the overall 
energy bias of the triple quantum dofP^. In this case, 
a “sweet spot” could be identified, where the noise is 
coupled only in second order to the qubitHowever, 
since each of the quantum dots is capacitively coupled to 
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FIG. 1. Three quantum dots in an (approximately) linear ar¬ 
rangement, coupled by virtual hopping between the quantum 
dots 1 and 2 and the quantum dots 2 and 3. Hopping be¬ 
tween the first and the third quantum dots is neglected due 
the nearly linear arrangement. Each quantum dot can be oc¬ 
cupied by a maximum of two electrons with opposite spins 
due the Pauli exclusion principle. 


its own electrodes, and to a different set of charge flnc- 
tnators in its immediate vicinity, the electrical potential 
on each quantum dot will in reality fluctuate indepen¬ 
dently, leading to a more general and more damaging 
noisy environment than previously assumed. The effec¬ 
tiveness of sweet spots has recently been investigated for 
a linear dot geometry in the special case of symmetric 
couplings ti = G, i.e., for y = O.^In this paper, we ad¬ 
dress the question whether coherent operation at a sweet 
spot of the RX qubit is still possible in these realistic 
conditions. To address this question, we study a noise 
model with each quantum dot coupled to independent 
stochastic charge fluctuations. 

This paper is organized as follows. In Sec. [TTj we in¬ 
troduce the model for the qubit under the influence of 
electrical noise investigated in this paper. Subsequently, 
in Sec. |III[ the noise model is derived and the dephasing 
times are calculated. We first show that a sweet spot, if 
it exists, needs to lie beyond the scope of a perturbative 
Schrieffer-Wolff (SW) approximation. To go beyond the 
SW approximation, we also include the low bias regime in 
our calculations. To this end, we calculate the spectrum 
of the exact Hamiltonian and find a sweet spot in a hy¬ 
bridized charge configuration of the triple dot. In Sec. |IV[ 
the pure dephasing times of the qubit at the sweet spot 
and in the RX regime are calculated in a Ramsey free 
decay setting. We conclude in Sec. |V| with a summary 
and an outlook. 


II. MODEL 


We consider three quantum dots with a single available 
orbital state in a linear arrangement (Fig. [^, described 


by the three-site extended Hubbard Hamiltonian, 


^Hub = ^ 


^ni{ni - 1 ) + ViTii 


+ 


E 

(bi> 


UcTliTlj + + h.C.^ 


a=n 


, ( 1 ) 


with the electron creation and annihilation operators cj ^ 
and Ci^a- on the site (dot) i with spin a, the number op¬ 
erator Ui = X]cr4crG,cr: thc pauwisc hopping matrix el¬ 
ements tij with i,j G {1,2,3} which can be controlled 
through variation of the gate voltages. Here, we con¬ 
sider symmetric nearest-neighbor hopping {tij = tji)^ 
and set tis = tsi = 0, t/ = V^ti 2 = t{l — y), and 
tr = \/2t23 = t(l + ^). Here, we have also introduced the 
mean hopping matrix element t and the hopping asym¬ 
metry parameter y. 

The first term in Eq. 0 describes the Coulomb energy 
U required when adding a second electron to one of the 
dots. The influence of the external gate electrodes is char¬ 
acterized by the second term Vifii, whereas Uc denotes 
the Coulomb interaction between electrons in nearest- 
neighbor dots. Filling the triple quantum dot with three 
electrons allows for 20 charge and spin configurations. 
There are two states each with a charge configuration 
(2,1,0), (2,0,1), (1,2,0), (1,0,2), (0,2,1), (0,1,2) and 
8 states with (1,1,1), where (m,n,/) denotes a charge 
state with m electrons on the left, n in the center, and 
I in the right dot. Defining the voltage difference 5 be¬ 
tween the outer dots as 5 = (Vi — V 3)/2 and the voltage 
difference Vm between the outer dots and the middle dot 
as Vm = (Vs + Vi — 2 V 2 )/ 2 , we note that charge transi¬ 
tions between (1,0,2) and (1,1,1) and between (2,0,1) and 
(1,1,1) occur at 5 = ±A = ±{U — 2Uc + Rm)- 

We work in a parameter regime of applied gate volt¬ 
ages Vi where only the charge states (1,1,1), (2,0,1) and 
(1,0,2) are accessible, and restrict ourselves to the sub¬ 
space with total spin S = 1/2 and spin 2 ; projection 
Sz = 1/2, spanned by the stated 




|0) — I'^)i3 1^)2 ~ ^ 

1^) 



ki,i/2) = k)n |t)3 = 4,t444,t 1^^^) ’ 

^3,1/2) = it)i 1^)33 ~ 4 ,t 4 ,t 4 ,t 1^^^) 5 (^) 


where |vac) denotes the vacuum state. Here, the states 
|0) and |1) are the logical qubit states of the exchange- 
only qubit in the (1,1,1) charge sector, while lsiq/ 2 ) 
and |si 1 / 2 ) denote the accessible states with the same 
spin but one doubly occupied quantum dot. This four¬ 
dimensional subspace can be separated from the remain¬ 
ing states by applying a large uniform external mag¬ 
netic field 5ext such that the states with spin projec¬ 
tion rus = ±1/2, ±3/2 along the ^ axis are split by the 
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FIG. 2. Schematic representation of a triple quantum dot 
confined by a potential V(x) occupied by three electrons (red 
dots). There are 20 possible conhgurations for the electrons 
if only one orbital per QD is considered. The hopping matrix 
element of electrons between the quantum dots is denoted by 
ti and tr. Here, we describe the response of the system to 
two noisy energy bias parameters, (a) the energy difference 
£ between the outer dots, with noise amplitude (5s, and (b) 
the energy difference A between the middle dot and the mean 
energy of the outer dots, with the noise amplitude (5A. 

Zeeman energy. The remaining = 1/2 states have 
either a total spin S' = 1/2 or S'= 3/2. The states with 
S = 3/2 and charge configuration (1,1,1) are almost com¬ 
pletely decoupled from the S = 1 /2 states if the exchange 
interaction is ongoing and much stronger than the Over- 
hauser field gradients.!^ Different states than the four 
defined above can be neglected if one assumes a strong 
Coulomb repulsion between electrons in neighboring dots 
(large Uc) and large energy gap between the orbital lev¬ 
els in such a manner that only the lowest orbitals are 
occupied. In the relevant subspace {0,1, 1 / 2 , <§ 3 , 1 / 2)5 
the Hamiltonian can be expressed as the 4x4 matrix 

/ 0 0 ti/2 tr/2 \ 

fr_ 0 0 V^ti/2 -V^trl2 

ti/2 V^ti/2 A + e 0 ■ 

\ tr/2 -^/Ztr|2 0 A-e / 

We use a simple model for the description of electrical 
noise in the gate voltages Vi by adding an independent 
noise term Wi to each gate voltage, = 4^ + Wi. Since 
only voltage differences affect the relevant spin dynam¬ 
ics, these three noise parameters can be reduced to two 
parameters. One noise parameter fe = {SVi — 5V^)/2 
represents noise in the voltage difference between the 


outer dots 5 , as shown in Fig. |^a). The second param¬ 
eter is (5A = { 6 Vi + SV 3)/2 — 6 V 2 and applies to the 
center voltage Vm = (V 3 + W ~ ‘ 2 y 2 )/ 2 ^ see Fig. ib). 
A third independent noise variable only leads to an in¬ 
consequential global phase of the qubit wavefunction. 
More precisely, we define the noisy voltage difference 
5 = {Vi — V 3 )/2 = eq 5e and the noisy effective center 
voltage A = U — 2Uc — Vm,o — = Aq — 6 A which 

also accounts for fluctuations in the Coulomb energies U 
and Uc- Here, the variables with subscript 0 indicate the 
parameters in the absence of noise. 

We define the RX qubit as the subspace spanned by 
the two eigenvectors corresponding to the smallest two 
eigenvalues of In its eigenbasis, the qubit Hamiltonian 
can then be expressed as 

^RX = 

with the Pauli cTz matrix and the energy splitting hw. 
The effect of the fluctuations 5e and 5 A in the parameters 
5 and A described by H in Eq. ^ leads to fluctuating 
terms in the RX qubit Hamiltonian. In the eigenbasis of 
the unperturbed RX qubit, these terms have the form 

Hrx = ^ ((^0 + 5^z)crz + SuJxCTx + Sujyay) (5) 

with the unperturbed eigenfrequency ujq and the longi¬ 
tudinal corrections (up to second order) 

5uj, = LOeSe + wa<5A + + uj^^aSsSA. 

(6) 

Here, the first derivatives of the qubit frequency uOp = 

I , determine the location of the sweet spot via the 
dp leoAo ^ 

condition = oja = 0 (see below), while the second 
derivatives ujp^q = with p^q = e, A limit the 

phase coherence of the RX qubit at the sweet spot. The 
transverse contributions SuJx,y are needed to calculate the 
qubit relaxation time, which is not our concern here. The 
longitudinal contribution 6 uJz represents the strength of 
the coupling between the qubit and the noise in first or¬ 
der and thus should be eliminated. Hence, the points 
(5, A) in parameter space where 6 uJz = 0 are known as 
’’sweet spots”. This becomes clear when expanding the 
eigenenergy difference from Eq. (§, 

U1 = (ulo + SuizP + Sul + 

— + 0{Slo^). (7) 

Away from the degeneracy lines 5 = ±A, the effect of 
the coupling to the (2,0,1), (1,0,2) states |si 1 / 2 ), 153 ^/ 2 ) 
can be taken into account using an effective Hamilto¬ 
nian in the low-energy (1,1,1) subspace which can be 
obtained by applying a Schrieffer-Wolff (SW) transfor¬ 
mation Hqr = He~^ such that the resulting matrix is 
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FIG. 3. Eigenfrequency Ct; as a function of either (a) £o with a minimum at £min ~ — (8/5)yAo or (b) A with a minimum 
at Amin ~ — (8/7)y£o. The parameter y — —0.1892 is derived from ti — 0.022 meV and tr = 0.015 meV^, while the other 
parameters are chosen to f = 0.018 meV, in (a) A = 3t and in (b) s = 3i. The Insets are magnifications of the area around the 
minima marked with a red rectangle. 


block-diagonal in lowest order S ^ ti^tr. As a result, we 
find the Heisenberg Hamiltonian, 


i^Heis = + (8) 


with the exchange energies Ji = / (A + 5 ) and Jr = 
tl/{A — e). In the logical subspace spanned by |0) and 
|1), the Heisenberg Hamiltonian becomes 

H = -^\0){0\-^- 11) (1| - (|0) {1| + |1) (0|) (9) 

with the mean exchange parameter J = {Ji + Jr)/2 and 
the exchange difference j = {Ji — Jr)l2. Diagonalizing 
we can write the RX qubit Hamiltonian in its eigenbasis, 
Eq. @ , with 


hiO = yfj‘^ + ‘dp. ( 10 ) 



EIG. 4. Three dimensional plot of energy gap 
The parameters are chosen to ti = 0.022 meV and tr = 
0.015 me V^. The clipped parts are diverging and outside 
the scope of the SW approximation. 


III. DEPHASING OF THE RX QUBIT 
A. Non-degenerate regime (SW approximation) 

In our analysis, we first investigate the special case 
of only one noisy detuning parameter, e.g. setting ei¬ 
ther SA = 0 or = 0. Our results in this simple case 
are plotted in Fig. and show a minimum of uj{s^ A) 
at (a) 5min ~ —(S/bj^Ao for fixed Aq and (b) Amin ~ 
— {S/7)yso for fixed 5o, where y = (tr —ti)/(tr de¬ 
notes the hopping asymmetry. These minima are sweet 


spots for one fluctuating parameter; one of them has been 
studied previousl}Ell. The qubit energy splitting is in gen¬ 
eral a function of both 5 and A, as shown in Fig. 

In the non-degenerate regime |A ± ^ we can 

study the influence of the electric charge noise on the RX 
qubit by expanding the low-energy Hamiltonian Eq. © 
to first order in Ss and ^A. We transform the result into 
the eigenbasis of the unperturbed Hamiltonian, leading 
us to Eq. (O, with 6 uJx = ^/3(J5j — j5J)l2uj^ 5ujy = 0, 
assuming that ti r are real-valued, and 6 uJz = —3( J6J + 
jSj)/2uj. Here, U = dJ/ds^.^AoSe ^ dJ/dA\,^^A^6A), 
and similarly for Sj. 
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FIG. 5. a) Plot of the energy gap between the lowest eigenenergies hw = E 2 — Ei as a function of the detuning parameters 
A. We set the hopping parameters to ti = 0.022 meV and tr = 0.015 meV. The black point marks the position of the sweet 
spot (0.0074 meV, 0.0015 meV), which is a saddle point. In b) the position of the sweet spot for the detuning parameters £ and 
A is plotted as a function of the hopping difference 5t — U — tr, with tr = 0.015 meV. 


In the high-bias regime, i.e., away from 5 = ±A, we 
can expect a sweet spot in the presence of both e and A 
noise if 6 uJz = 0, i.e., if the derivatives of ujz with respect 
to both parameters A and 5 vanish. We find 


duj If ^dJ ^ .dj\ 1 1 


de (jo \ de de J uj /S? — e‘^ 


(ew^ - 4JjA) , 

( 11 ) 


dw f dJ o • W 1 1 
dA ujK dA ^dAj uj A^ - 


(-Aw^ + AJje) , 


where we have used 

dJ dj 1 1 


de dA w A2 - e2 
dJ dj 1 1 


(eJ-Aj), 


dA de w A2 — e 2 


i-AJ + ej). 


( 12 ) 

(13) 

(14) 


The condition ^ = 0 cannot be fulfilled for 5 7 ^ 

A. Therefore, we cannot find a sweet spot in the non¬ 
degenerate regime. 


B. Degenerate regime (exact solution) 

Since no sweet spot is found within the scope of the SW 
approximation, we now investigate the degenerate regime 
|5 ± A| ^ti^ri which is outside the scope of the SW ap¬ 
proximation, and in particular the low bias regime A, 5 ^ 
0. For this purpose, we directly calculate the eigenener¬ 
gies of the subspace Hubbard Hamiltonian, Eq. (§. We 
denote the eigenenergies Ei < E 2 < < E 4 and note 

that they are functions of the two detuning parameters 5 
and A as well as the two hopping parameters ti and G. 
Analytical expressions for Ei can be obtained, and are 


shown in Appendix The qubit states are defined as 
the two lowest energy levels which match the RX qubits 
states in the (1,1,1) charge sector, with energy separation 
fuj = E 2 — El. In Fig. I^a) we plot huj for fixed hopping 
parameters. Here, we indeed find a sweet spot (indicated 
with a black dot) near but not exactly at 5 = A = 0. 
The position of the sweet spot in ( 5 , A) space is shown 
in Fig. I^b) for G = 0.015 meV as a function of the hop¬ 
ping strength difference 6 t = ti — E. The formula used 
to calculate the energy gap and the sweet spot is given 
in Appendix The resulting sweet spots always fulfill 
the condition 5 > |A|, hence they are located outside the 
strict (1,1,1) charge configuration and the qubit states 
acquire a component of states with a double occupation 
of the right dot (1,0,2) and the left dot (2,0,1). However, 
being a sweet spot, the qubit at this working point is 
only weakly coupled to charge noise. In the special case 
of symmetric hopping, ti = tr, we find a sweet spot at 
e = A = 0 . One could expect that leakage is problem 
in this case since the energy gap between the two logi¬ 
cal qubit states is comparable to the energy difference to 
other non-qubit states, whereas in the RX regime the en¬ 
ergy gap between the logical qubits is far away from other 
states. However, the dynamics show greatly suppressed 
leakage if only one parameter 5 or A is driven with the 
resonant frequency ci;(5,A). E.g., Rabi transitions be¬ 
tween the qubit states are much faster if 5 is periodically 
driven in contrast to transitions between the energy lev¬ 
els E 2 , E 3 , where the detuning parameter A needs to be 
driven. The reason for the sensitivity to only one pa¬ 
rameter is the symmetry of the energy difference with 
respect to the detuning parameters. In the fully sym¬ 
metric case ti = tr, Rabi transitions between the energy 
gaps only occur for one driving parameter, while driving 
with the other parameter is completely suppressed. Eor 
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asymmetric hopping the symmetry is weakly broken and 
Rabi transitions can occur for both parameters but with 
varying speed. 


IV. PURE DEPHASING 

In this section we investigate the effects of charge noise 
on the RX qubit at the sweet spot in the asymmet¬ 
ric charge configuration (degenerat e regi me) in compari¬ 
son to the effects on the RX-qubits^^^ within the non¬ 
degenerate (standard RX regime) with symmetric (1,14) 
charge configuration. Since no sweet spot can be found 
within the RX regime the eigenenergies couple linearly 
to both noisy parameters Ss and SA giving rise to a de¬ 
phasing times which scales inversely proportional with 
the noise amplitude. However within the RX regime, one 
can find the best working points, where one parameters is 
minimized, e.g. A ^ —^jlye or 5 ~ —(8/5)^ A, which 
corresponds to a sweet spot for one parameter, where the 
dephasing time is inversely proportional in the other pa¬ 
rameter. At the real sweet spot found in this work this 
scaling is at least inversely quadratic. This characteristic 
trait can be observed in Fig. (a), where the estimated 
dephasing time is plotted as a function of the n oise l evel. 
For a current noise level on the order of the 

RX regime appears to be the better choice, since the re¬ 
sulting dephasing times in the RX regime are two orders 
of magnitudes longer than at the sweet spots. However, 
below a noise level on the order of 0.1 /reV, it becomes 
advantageous to choose the sweet spots due to the better 
scaling. In the subsection below, we resent the free decay 
model used for the calculation of T^. 


A. Dephasing model 

To study dephasing, we start from the noisy RX qubit 
Hamiltonian Eq. (§ and focus on the longitudinal noise 
SuJz. The time evolution operator U{t^tQ) from an initial 
time to to some later time t can be written as 

U{t, to) = exp [-i(p{t)az] (15) 






1 

71 


|0)+e'^W|l) 


(17) 


One observable of interest is the mapping on the initial 
state P = |-7) (+1 which leads to the free decay ansatJ^ 


{P) = \{+\U{t,to)\+)f 


! + /(*) 


(18) 


The function f{t) describes the dephasing in a free decay 
model and is given for Gaussian distributed noise by 


f{t) = (e‘'^> exp 



(19) 


The detailed formula for the decay and the derivation can 
be found in Appendix [B| For t ^ oo the superposition 
is destroyed and the expectation value is 1/2 for both 
states, as expected. 


B. Approaching real systems 


For further calculations, such as evaluating the inte¬ 
gral in Eq. ( p!^ , we require the knowledge of the power 
spectral density S{(jj) of the noise, hence we have to 
consider electric charge noise in a more detailed man¬ 
ner. Here, we consider Gaussian distributed noise with 
a power spectral density S{(jj) = A\uj\~^ with variance 
A = of the noise Se (SA) and 7 = 1 which resem¬ 

bles charge noise in double quantum dots.l^ The analysis 
of Eq. ( [T^ leads to Gaussian behavioil^ for the decay 

rate, fi{t) oc exp — 0{t^) , with 


=h 


Y Ae logr + ^ Aa logr 




+ (1 + 2 logr) + 

2 ^ -j- UJs^sUJa,A 




As A A 


Aa (1 + 2logr) 

( 20 ) 


with the accumulated phase 


0(t) = f dt'SuJz{t') 

Jtn 


[ dt' 

Jto 


UJsds(t') + UJA^A(t') + —UJs,£d£(t')‘^ 




(16) 


The time ordering operator is not needed, because only 
longitudinal coupling {az) is considered. Considering the 
effects of an initial coherent superposition of the qubit 


with r as the quotient of upper and lower cutoff, which 
is needed for ensuring convergence of the integral. The 
detailed derivation of the formula above can be found in 
Appendix B. Since at the sweet spot the terms linearly 
coupled to the noise vanish, one should expect long de¬ 
phasing times T^. However, our findings shown in Eig. 
(e) exhibit shorter dephasing times in comparison with 
the RX regime (Eig. |^(c), area with long times) due 
to a strong contribution of the second order couplings. 
This contribution, especially ujs,£^ strongly limits the de- 
phasin g tim e at a realistic noise level in the order 
of /ieV®^. Small improvements can be made by con¬ 
sidering larger hopping parameters because ujs^s ^ ^ 
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FIG. 6. (a) Doubly logarithmic plot of the dephasing time as a function of the rms (root mean square) value of the noise 
at the best working points (Sec. Ill A) in the RX regime (red) and at the sweet spot (blue). The parameters are chosen as 
A = 0.32 meV, ti = 22 /xeV, U = 15 /ieV and e: = — (8/5)yA for the RX regime plot and £ = 74 /ieV, A = 15 /xeV, ti = 0.22 meV 
and tr = 0.15 meV for the plot of the sweet spot, (b)-(e) Density plots of dephasing time as a function of parameter space 
for different noise amplitudes (black circles in (a)). Note that for the lower plots the hopping parameters are ten times greater 
than for the upper plots. The black points indicate the position of the best working points plotted in the left (upper) and the 
sweet spot (lower). The black squares in the upper plots marks the space plotted in the lower ones. 


or (to a small extent) with a stronger asymmetry, but 
nonetheless, the dephasing times remain several orders 
of magnitude shorter than the best points of operation 
within the RX regime. Overall, for the currently avail¬ 
able noise level, we find that the sweet spots do not lead 
to an improvement in coherence and the standard RX 
regime should be favored instead. Those best operation 
points can be achieved by fine-tuning 5 and A in such a 
manner that either of the two parameters is minimized 
(typically 5 ), while staying within the (1,1,1) charge con¬ 
figuration regime. This limits the pure dephasing time 
to a maximum at 5 ~ ±4 meV since overly large param¬ 
eters A and 5 give rise to other charge configurations, 
effectively reducing the benefit gained by leaving the RX 
regime.^ Importantly, the situation changes completely 
when lower noise levels become available, because of the 
quadratic scaling behavior of the dephasing times at 
the sweet spot compared to the linear scaling of in the 
RX regime. This different scaling is outlined in Fig.[^(a). 
For a noise level of 10“^ /icV we find at the sweet spot 
one order of magnitude greater than in the RX regime 
(Fig. § (b) and (d)). The crossover between the two 
regimes occurs at approximately one order of magnitude 
less than the currently measured noise leevels. Hence, 
by purifying the materials or improving the noise filters 
in such a manner that the charge noise level is lowered, 
these sweet spots offer a promising perspective in further 
reducing dephasing of charge noise in the future. 


V. CONCLUSION AND OUTLOOK 

In this work, we have presented a full description for 
charge noise in the RX qubit. We have shown that there 
are sweet spots for two different coupled noise parame¬ 
ters, which are suppressing charge noise coupled to the 
RX qubit in first order and give rise to operation points 
with quadratic noise terms. However, considering both 
noise parameters simultaneously, no suitable sweet spot 
is found within the scope of the SW transformation re¬ 
sulting in qubit states linearly coupled to noise. By tak¬ 
ing into account the low bias regime, we found a sweet 
spot outside the scope of the SW approximation in the 
crossover region to the (2,0,1) and (1,0,2) charge configu¬ 
ration, with a precisely determined location in the ( 5 , A) 
parameter plane depending on the hopping asymmetry. 

For the description of the dephasing of the RX qubit 
we used a Ramsey free decay model to describe the re¬ 
sulting dephasing times T^. We also included quadratic 
effects which dominate the dephasing at the sweet spots. 
As a result, we found a Gaussian behavior of the dephas¬ 
ing as a function of time in lowest order. In the next step 
we compared the usefulness of the sweet spots with the 
best working points within the RX regime. The result¬ 
ing analysis shows that the best working points within 
the RX regime should be favored for currently available 
noise levels. However, if we consider an improvement by 
about one order of magnitude in the charge noise level. 























our sweet spots should be the favored option due to their 
better scaling behavior of the qubit coherence time. In 
this work we have neglected the influence of other noise 
than charg e noise such as spin orbit coupling, hyperfine 
interaction,fluctuations of the homogeneous mag¬ 
netic field and so on. In future studies, these effects can 
be included in a full quantum master-equation approach 
to further improve the results. 
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a =A2 -e^ - tf - tf. (A2) 


In our investigation we only consider the energy gap be¬ 
tween the lowest two levels u = Ei — E 2 = k,- and its 
partial derivatives and The condition for the 

sweet spot = 0 is giving rise to a system of 

nonlinear equations. Instead we minimized the equiva¬ 
lent system 


Appendix A: Exact solution of the subspace 
Hubbard Hamiltonian 



(A3) 


To describe the parameter space for small 5, A, which 
is outside the scope of the SW approximation, we calcu¬ 
late the eigenenergies of the Hubbard Hamiltonian in the 
subspace spanned by {|0), |1), |si,i/ 2 ), IS 3 . 1 / 2 )}, Eq. 
directly with the general solution for a polynomial of or¬ 
der four. This is giving rise to the four eigenenergies 
(numbered from lowest to highest) 


using that all parameters are real valued. 


Appendix B: Derivation of the free decay rates 


For an estimation of the dephasing rate, the observable 
of interest is the projection on the initial state Eq. (18). 
Therefore, it is sufficient to calculate 


El 

E 2 


Es 


E 4 



m = (e-q, (Bi) 

which can be expanded by using the cumulants 

log(</.(/))“ + •••. (B2) 

Assuming Gaussian noise with zero mean {Ss{t)) = 
{6A{t)) = 0, all cumulants higher than two and all odd 
moments vanish, resulting in 
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log(e“^) = - i Us5e{t')dt^ ) "^ ( (/ ^ ^e,e5e{t'fdt' 


+ 4 




+ U)e,ede{t')'^dt^ (^j LOA,AdA{t')‘^dt' 


Fourier transforming and some calculus leads 


log(e“^) = - j dA Se{uj)smc^ - 


ujt\ uj\ 


I sine 


^ i -2 
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duo duo'sinc^ 


duo duo'sinc^ 
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uo -\- uo' 
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duo Sa{^) si 

-oo 

5,(a;)5,(a;')) - 
5A(a;)5A(w') 


u:s,eU:A,At {SS^){SA^) 

-^-k^{Se^}{SA^), 


(B3) 


(B4) 


where Sq{uo) = f^^Sq{r)e~^^'^dr {q G { 5 , A}) denotes 
the spectral energy density of the noise coupling to bias 
parameter q. 

1. First Order 


which diverges at the lower limit. To ensure convergence 
of the integral the spectral density is modified to cutoff 
the lowest {uo < uor) and highest frequencies {ub > uou) 
with uou > 00 R > 0, hence S{uo) = j^Q{ub — uor)Q{uou — 
ub). This modification leads to 


For our investigation we consider noise with a spectral 
which is anti-proportional to the 


density S'(w) = 
frequency. Set into equation (B4) yields the free decay 


f{t) = exp 
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dd) Tginc^ (— 
A V 2 


(B5) 


fi{t) = exp 


2 2 . r’-' ^ . 2 /At 

-2t ojgAg J dw y sm (^y 


which can be integrated to 


(B6) 


flit) = exp 


-2Wq Ag 


+ - 


— 1 + cos{tuJii) + t^uj^ciitujfi) — tujfismitbjfi) 


4^1 


1 + cos{tuou) + t‘^uo^ci{tuou) — tuou sm{tuou) \ 

4w/ ) 


(B7) 


with cosine integral ci(x) = — dt. In the limit uor oo, all high frequencies are valued, the second term in 

the exponent vanishes and the solution is 


7 \ r, 2 A f -1cos(twfl) + - twflsin(twfl)y 

JP,oo{t) — exp ZiOqAq ^ ^^2 j 

Considering small tuoR 1, the first order of the evolution of the exponent can be written as 


/i,oo(i) ~ exp 


OJ^ A. 
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log 
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Since only the lower frequency regime is dominated by pink noise and pink noise can be neglected for higher frequencies, 
a high frequency cutoff can be motivated and the two terms in equation (B7) can be simplified to 


1 / 1 


1 


t 




/i(i) = exp -2ujgAq( - { — - 2 ]+T t) - ci(w_R t)) + — (sinc(wfl t) - smc{uju t)) 

\ 4 V LO^ 


=const>0 

Considering both tojR 1 and toju :§> 1 to be small, the exponent can be expanded and rewritten a^^ 


(BIO) 


fi{t) « exp 




(X exp 


2/1^1 f 

^ 4 \ujR 


(Bll) 


Introducing r = ^ the expression above takes the form of the first two terms in Eq. (20). 


2. Second Order 


For the second order terms we have to calculate 
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duj duo'sinc^ 


UJ uj' 


Ss,a{^)Ssa{^')^ 


(B12) 


from Eq. ( |B4[ ) which can be done in a similar manner than the linear term above by introducing the same power 
spectral density S{uj) = ^©(o; — ujr)Q{uju — dj) with a low and high frequency cutoff. The integration gives 
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but can be approximated by a Taylor series expansion into the familiar expression 


(B13) 
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With ^ ^ the expression above takes the form of the missing two terms in Eq. (20). 
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